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Houdayer and Martin reply: 

Marinari, Parisi, and Zuliani have studied the 
Edwards- Anderson spin glass model at a field B = 0.4. 
They had previously studied the four-dimensional case ffl] 
and in their comment [Q to our paper || they present 
results for the d = 3 case using the same techniques. The 
cornerstone of their approach is an out-of-equilibrium es- 
timate of <7 mm . They find this quantity to be different 
from the (equilibrium) mean value of q, giving evidence 
for replica symmetry breaking (RSB). However we see 
a danger in relying on out-of-equilibrium measurements: 
metastable states that do not contribute to the (equilib- 
rium) P(q) (because they have excess free energies di- 
verging with system size) may very well contribute to 
out-of-equilibrium overlaps. 

Since we have some doubts about the validity of out- 
of-equilibrium measurements of 9 mm , let us consider the 
evidence for RSB in the presence of a field using equilib- 
rium measurements. Most work has been performed in 
d = 4 where there is a clear signal of a growing spin glass 
susceptibility. Less clear is whether this quantity actu- 
ally diverges at T > and B > 0: fits are compatible 
with such a divergence but it is difficult to conclude that 
there is a finite temperature transition. Measurements 
of higher order cumulants of P(q) are very disappoint- 
ing: the cumulants do not cross as they do in zero field, 
there is no clear pattern in the data, and in fact there is 
no sensible way to extrapolate the data to larger sizes. 
Furthermore, P(q) has a long tail at q < that cannot 
be there in the thermodynamic limit, and there is no hint 
yet of a delta function peak at q = <7 mm - The situation in 
d = 3 is even more ambiguous because the evidence for a 
diverging spin glass susceptibility is much weaker. Never- 
theless, since simulations of the Sherrington-Kirkpatrick 
model run into similar difficulties, one need not conclude 
that such results disfavor a RSB scenario in finite di- 
mensions. But it is also fair to say that there is today no 
substantial evidence via equilibrium measurements for an 
Almeida-Thouless transition line in d = 3. 

Part of the difficulty stems from the nearby critical 
point (T c , B = 0) that leads to severe finite size effects; 
this may explain the non crossing of the Binder cumu- 
lants for different size lattices. Staying away from that 
critical point requires lowering the temperature, leading 
to insurmountable difficulties for thermalizing the lat- 
tices. Our approach bypasses this problem by taking the 
zero temperature limit and finding ground states. Do- 
ing so, we found that the finite size effects for the mean 
field model were very small as shown in figure 3 of our 
paper. A comparison to what occurs in simulations at 
finite temperature of the Sherrington-Kirkpatrick model 
suggests that the (T c , B — 0) critical point is quite far 
away from where we work. Extrapolating this to the 
Edwards-Anderson model case, we are led to conclude 
that the crossing points of the curves of figure 1 in our 
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FIG. 1. r(N,B,0.15)/VN versus B. 

paper simply converge to B = B c = 0. We tried to sub- 
stantiate this hypothesis by finite size scaling (figure 2 in 
our paper). If on the other hand one insists on having 
a critical value B c > of the field, we are led to ask 
whether the curves r(N, B, 0.15) /y/N vs B superpose at 
large TV on a curve extending to B > 0. Our data is dis- 
played in figure |l|. Judging from this figure, we probably 
would need to work with lattices larger than 20 3 before 
such behavior could be seen. This may be feasible in the 
not so distant future, but much remains to be done. 

J. Houdayer and O.C. Martin 

LPTMS, Universite Paris-Sud, F-91405 Orsay, France 

July 22, 1999 

PACS numbers: 75.10.Nr, 64.60.Cn 
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